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Abstract 

We study SU(iV) level k Chern-Simons theories coupling to a fundamental 
fermion with chemical potential for U(l) flavor symmetry. We solve this system 
exactly in the 't Hooft limit, N, k — > oo with A = N/k fixed. The solution exists 
up to |A| = 1 for a fixed value of chemical potential. We study the system 
beyond |A| = 1 by considering U(iV) gauge group and taking the diagonal U(l) 
holonomy into account. The result suggests the system becomes unstable and 
a phase transition happens at |A| = 1. 



1 Introduction 



Vector (or sigma) models play important roles to understand not only phenomenology but 
also theoretical high energy physics. In the context of AdS / CFT duality they play a crucial 
role for the duality proposed in The duality claims that the singlet sector of a three 
dimensional O(iV) (or U(N)) vector model with conformal symmetry and a higher spin 
gravity theory on AdS4 J2|-|J in the large N limit. This conjecture has been tested by 
comparing spectra and correlation functions of both sides |5|-f§ . 

Motivated by this conjecture three dimensional vector models have been studied in 
presence of CS interaction by gauging global symmetry |9|, [T0|] . It was studied that confor- 
mal symmetry and higher spin symmetry are still preserved up to 1/N correction and the 
effect of CS interaction is reflected as a suitable change of boundary condition for a bulk 
scalar field. This class of conformal field theories turned out to have strongly constrained 



correlation functions \~. 4|. The duality can be also studied by changing a base manifold 



15|, |l6ll or in accordance with supersymmetry || |T7Hl9H • A recent review is |20| 



In this paper we consider a little generalization of by adding chemical potential for 
U(l) flavor symmetry into Chern-Simons-fermion vector model. One of the achievements 
in || is solving the system by taking the light-cone gauge and the 't Hooft limit. The 
free energy is determined exactly as a function of A. Interestingly the solution exists up 
to |A| = 1 and the free energy vanishes at that point. The motivation of this paper is 
to compute other various thermodynamic quantities such as the grand canonical potential 
and the thermodynamic descendants thereof exactly as a function of A and to address what 
happens to the system beyond |A| = 1. 

The rest of the paper is organized in the following manner. In Section |^ we present 
the exact solution of Chern-Simons-fermion vector model with chemical potential in the 
't Hooft limit. In Section ||| we calculate various thermodynamic quantities exactly as a 
function of A. In Section | we study this system beyond |A| = 1 by considering the gauge 
group U(iV) and taking account of the diagonal U(l) holonomy. The final section is devoted 
to summary and discussion. 



2 Chern-Simons-fermion vector theory with chemical 
potential 

In this section we study SU(iV) Chern-Simons gauge theory at level k with a fermion which 
is in the fundamental representation of the gauge group in presence of chemical potential 
fi for U(l) flavor symmetry. The action of this system is 



S 



d 3 x 



(2.1; 



We use the same convention in [19 



This sytem can be exactly solved in the 't Hooft limit || 

N, k — > oo with A := N/k fixed. 



(2.2) 
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The method is the following. Firstly one fixes the gauge degrees of freedom by the light-cone 
gauge in Euclidean space A_ := ^75(^1 + 1A2) = and integrate out the other components 
of the gauge field. This generates four fermi interaction. Secondly one introduces a bilocal 
singlet field S so that the four fermi interaction term cancels, which is known as Hubbard- 
Stratonovich transformation. S turns out to coincide with the self energy of the fermionic 
field on shell. The resulting action is quadratic with respect to the fermion and thus we can 
integrate it out by gaussian integration. The remaining action is written in terms of the 
singlet field and N can be factored out from the action. Finally one can obtain the exact 
effective action in the leading of 1/N expansion by the steepest descent method under the 
assumption of translational invariance. 

From the action fl2.1|) we can easily see that the result in our case can be obtained from 
that in || by exchanging the momentum in the fermion propagator in such a way that 

P^-^P^-=P^ + i^^ 3 . (2.3) 

After this repalacement we find the exact effective action 



log + £(<?)] + ^£(g) 



(2.4) 



2 V*7^m + s (<?), 

where the trace is taken over a two-by-two matrix with spinor indices and V is the volume 
of the Euclidean space. As shown in || this effective action correctly reproduces the 
connected vacuum diagram consisting of planar graphs perturbatively with respect to A. 
The exact propagator of fermionic field is 

(<K-2#(g)> = . ~\ T (2n) 3 5 3 (-p + q). (2.5) 

Spinor indices are suppressed here and hereafter. We also find the saddle point equation 

^ f-.+ 1 .3 3 1 „+\ 27TA 



E ^ = / 7Ti3 7 + . ^ s t - r . M ~ , w , r - ( — r . (2-6) 

J (2tt) 3 V 27%* + «7 M 9/i + s (s) /KP-I)- 

where p± = 775 (f>i T 7 = 771(71 ^ 272)- We can check that this equation precisely 
agrees with Schwinger- Dyson equation in the 't Hooft limit and also coincides with the 
equation of 1PI self-energy diagram, which is now given only by the rainbow diagram ||. 
Let us solve the saddle point equation ( |2lf ). For this end we decompose S into 

£ = 7 "E„ + IS/ (2.7) 

where / is a two by two unit matrix. Plugging this into (|2.6|) we find 



By comparing both sides we obtain 

S_ = S 3 = 0, (2.9) 

4ttA /• A / zg_ \ 1 

1{p> i J (2^\(q,-z^(q)y + j: 1 (q) 2 J(p-q)- [ ' 1 

y M = / _gg_ ( ^M) \ 1 (2U) 

+ {P) i J (27r)3V(g„-iE„(?)) a + Si(g)V(p-?)- l ' J 
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From the dimensional analysis and rotational covariance in two plane we can set 

S/(p) = f(p')p s , S + (p) = ig(p')p + (2.12) 

where / and g are unknown functions of p' = — . Plugging these into the above equations 
we obtain 



f(p') Ps = 4vrA J 

g{p')p+ -- 



4ttA 



d 3 q 

V (i + g(q'M + oi + (f(q')Qs) 2 J (p - g)_ 



+ <?(?'))<7 S 2 + 51 + (/Ma.) V (p - ^) 



(2.13) 
(2.14) 



By using 
to 



_9 1_ 

9p+ (p-<?)- 



27r5 2 (p — g), we can show ((f(p')p s ) 2 + g(p')p 2 ) = 0, which leads 



c 



f(p'Y + g(p') = ^ 



p 



(2.15) 



where Co is a dimensionless constant which will depend on A. g can be determined by ( p. 15 ) 
so we focus on solving the equation of /, ( |2.13p . Substituting ( |2.15|) into fl2.13| ) we obtain 



f(p')p s = 4ttA J 



d 3 q 



Q- 



(2tt) 3 \q 2 s +qi + c fi 2 J (p - q). 
Integrating the angular part of p_ we get 

q s dq s f dq 3 -1 



f(p')Ps = 4vrA 



27r J 2tt q 2 s + gf + c /i 



(2.16) 



(2.17) 



We regularize UV divergence by dimensional regularization, as done in || |19[], which leads 

toQ 



f(p') = ^V\p' 2 + Co\0(V\p' 2 + Co\ 



1) -9(1 



\p ,2 + c \ 



(2.19) 



where 9(x) is a step function taking unit value when x > otherwise zero. 

Co can be determined as follows. First let us plug ( f2.19| ) into ( |2.15| ) and expand g 
around p' ~ 0: 



A 2 (vW( - 1) - 9(1 - v^ol) 



P 



12 



(2.20) 



From ( |2.14 ) and ( [2.19 ), on the other hand, we can show g(p') ~ around p' ~ in the 
same way as in |9|, [l9j . Therefore Cq has to satisfy 



Co - A 2 (v^o>(v^oT " 1) " 0(1 " V^oT)' 



■""We note a usuful formula to derive this equation 

dq3 



r (93 + «/-i) 2 + a 2 \a\ 



-0(|a|-|/i|). 



(2.21) 



(2.18) 
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A unique solution exists only when |A| < 1 as Co = A 2 and no solution when |A| > 1. 
As a result we obtain the solution of the saddle point equation 

f(p') = ^VT+A^^VT + ^-l) " 0(1 " a/p^TP)), (2.22) 

g{p') = ^(1 - p' 2 - A 2 MV7?P - 1), (2.23) 

where |A| < 1. X finally becomes 

E(p) = f(p') Ps I + ig{p')i + p+. (2.24) 

Remark that the effective mass becomes |A/x|. 

On the solution ( p.24|) the effective action (|2.4|) can be evaluated by using the dimen- 
sional regularization as 

Seff/V = -iV A2(1 ^ |A|) H 3 . (2.25) 
We can also evaluate momentum distribution of U(l) flavor charge density 

(nps) := SW?S W) (^ (p)7 ^ (g) > - (^Ph'Ml)) U=o) , (2.26) 

where the normalization is such that it vanishes when the chemical potential does. By 
using (|2.5|) and (|2.24|) this becomes 



(n Ps ) =N(l- 6{^vA + JM 2 ~ I*) ~ 6(-Vp 2 s + M 2 ~ A*)) • (2-27) 

This momentum distribution indicates this system is described by an ideal Fermi gas whose 
Fermi surface is p s = — A 2 . 

To summarize, the Chern-Simons-fermion vector model in presence of chemical poten- 
tial fi is effectively described by a free fermion system with the mass squared (A//) 2 under 
the 't Hooft limit .0 



3 Finite temperature analysis 

In this section we study the Chern-Simons-fermion system in nonzero temperature T by 
placing the system on R 2 x S 1 with circumfirence /3 = 1/T. We impose anti-periodic 
boundary condition for S 1 direction on the fermionic field in order to study the (grand) 
canonical ensemble thereof. Since Fourier mode in the thermal direction is quantized by 
P3 = 27r ^ + 2 - ; we can obtain the result in nonzero temperature from the previous calculation 
by doing replacement such that 

J H pz-F H neZ ' 

2 In this specific model we analytically see an important assumption in statistical physics that a weakly 
interacting many body system reduces to a free ideal one with appropriate parametrization in the ther- 
modynamic limit pl|]. 
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where F(ps) is an arbitrary function of the momentum in the thermal direction. 
We first solve the Schwinger- Dyson equation, which is now given by 



1 q n 1 , \ 27rA 

7 + - „~ , w x 7 3 " 7 3 . M ~ ^ w , 7 + T7 r • 3.2 

ij^ + E (g) z7^g M + E (g) Ji[p-q)- 



Repeating the same calculation we obtain 

S_ = E 3 = 0, (3.3) 

= 4vriVl v / rf 2 g / zg- \ 1 , 

^ ^£^7 (27r) 2 V(g-iE(g)) 2 + E 1 (g)V(p-g)_ 

4W1 v y rf 2 g / Ei(g) \ 1 

«'* (2 7 r) 2 V(g-^(g)) 2 + E 1 (g) 2 y'(p-g)_ 1 " ' 

From the dimensional ground the ansatz is such that 

^i(p) = f(Pifr)Ps, E+(p) = ig(p,jl)p + (3.6) 

where / and g are unknown functions of p = jr and jl — ^. Hereafter we will suppress /t 
in the argument. Then the above equations become 

q- \ i 



fip)Ps S ^ ( (1 + 9mq " + ^ + C^k)V (p - ?)- 

g(pK = -4^1 gm J 8 +£ + mqa) *) ( 3 - 7 ) 



The same argument in the previous section leads to 



c 



/(p) 2 + <?(p) = ^, (3.* 



where c is a dimensionless constant depending on A and /t. Substituting fl3.8| ) into (|3] 
we obtain 

/(p) P , =4*Ai £ / ^ ( g , + |- +cT2 ) 



10 g.dg s 1 ^ -1 



=W ^E^n^- (-) 



Regularizing UV divergence by dimensional regularization we findQ 

A 

V 



f(p) = b. log (^2(cosh ^\p 2 + c| + cosh/i)) . (3.11) 



3 A useful formula to derive this is 

taml ^) + tanh£M 



1 y 1 = tanh^™+tanh^y^ 
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c can be determined in the same way as in the previous section. Substituting ( p. lip 
into (|3.8| ) and expanding g around p ~ we find 



sip) ~ ^ 



A 2 ^log ^2 (cosh \J\c\ + cosh fx) 
p2 



(3.12) 



We can show g(p) ~ around p' ~ from (|3.7|) and ( p.ll| ). Therefore Q3.12j ) implies c > 
and 



a/c = |A| log (2 (cosh yfc + cosh /a)) 



(3.13) 



A solution exists when |A| < 1 and does not when |A| > 1. This fact can be easily seen by 
rewriting ( |3.13|) as 

(3.14) 



which suggests ei A i > and thus |A| < 1. 

Let us summarize the solution of the Schwinger- Dyson equation 

fip) = 4 log (2(cosh a/p 2 + c + cosh/2)) , g(p) = - /(p) 2 , 
p \ / p z 



(3.15) 



where a/c is determined by ( |3.13| ) or ( |3.14| ). Note that a/c physically means the thermal 
mass of the fermion divided by T. 

As in the previous section we can compute the momentum distribution of U(l) charge 
density in finite temperature 



(3.16) 



where we did replacement ( p.l|) in (|2.26|) . Substituting the above solution and simplifying 
it we find 



(rh.) = N - 



tanh vt±L±. + tanh 



-y/p 2 +c+A 
2 



(3.17) 

This is the momentum distribution of the free fermion gas in nonzero temperature with 
thermal mass a/cT. The fermi surface clearly seen in ( |2.27|) is now smeared by finite 
temperature effect. 

Let us compute the grand canonical potential density in this system, which is defined 

by 



G(T,li) 



S, 



eff 



s, 



e//|7>=0 



V 2 (3 



(3.18) 



where V2 is the volume of two plane given by V-i = V/ (3 and S e ff is the effective action on 
R 2 x S 1 with chemical potential 



Seff = NV 2 

qs-F 



d 2 q 

(27T) 2 



tr 



log + + 



i^q^ + S(g) 



(3.19) 




and we normalize the grand canonical potential as zero when T = ji = 0. By employing 
the dimensional regularization the grand canonical potential density can be computed asQ 

]\Tn~>3 / /~3 roc poo \ 

G(T,fi) = — — ^--^ + 3 J dyy]n(l + e-*-»)+3 J ^ dy y\n (l + e~^)\ . 

(3.21) 

We plot y/c and — G |,^ as a function of A varying fi in Fig. |l(a)| , Fig. [l(b)| , respectively We 
have assumed A > here. Those figures indicate that as |A| goes larger the thermal mass 
becomes bigger and the gratitude of the grand canonical potential smaller, and finally when 
| A| reaches one the thermal mass blows up and the O(N) grand canonical potential vanishes. 



Note that this grand canonical potential ( |3.21| ) reduces to the free energy density obtained 
in when fi — > and to the effective action ( |2.25|) when T — > since y/c — > |A/x|/T from 

The behaviors of y/c and G(T,fi) around p, ~ oo, which means low temperature or 
small chemical potential, with A fixed is the following. 



Xp, + \e {x - 1)fl + ■■■ , (3.22) 

6tU (1 - aV+ t 



NT 3 f\ 7T 2 

G(T, n) — -(1 - A 2 )/i 3 + — p, + • • • ) . (3.23) 
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We note a useful formula 



1 £ l0 S (( Mn o ^ + + h 2 )= log(l + e~^ h+ ^) + log(l + e -^ fe -M)) + ft. ( 3 .20) 



^ riGZ 
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(a) 



(b) 



(c) 



Figure 2: The U(l) flavor charge density, entropy density and heat capacitance normalized 
by NT 2 are plotted clS cL function of A varying fi/T = 0,1,2,3,4 from the bottom in 
Fig. |2(a)| , Fig. |2(b)| , Fig. p(c)| , respectively. 



Their behaviors around A ~ 1 with ft fixed are 



,'2cosh/t\ /2cosh/T 
V c ~ log — r- - log log I — — 1 - 



1-A 



1- A 



, iVT 3 (l - A) / r3 3 r 2 



Qtt 



cosh ft 



(3.24) 
(3.25) 



From the grand canonical potential one can obtain its thermodynamical descendants 
such as the pressure p, the mean U(l) flavor charge density Q, the average internal energy 
U, the entropy density S and heat capacitance C with fixed two dimensional volume. 

p = -G,Q = ~(d,G) T , U = dpifiG)^ S=- i^—j ,C=-T l^j , (3.26) 

where subscripts are such that partial derivatives are taken with them fixed. From a 
straight forward calculation we obtain 

NT 2 NT 2 a NT 2 2 u 2 

Q = -^-9, U = -2G, S = ^-{g-^gf), C = —(g - \ftg> + ^-g"), (3.27) 

where we set g = — j^G{T, //) = g(ft) and g' = = ^f-Q In particular we can show 

U and Q are related in large ft region by U ~ w gN }il\2\ Q^ ■ 



4 Comment on instability in |A| > 1 



In the previous section we solved the Chern-Simons-fermion vector theory with chemical 
potential in finite temperature in the 't Hooft limit. The solution exists when |A| < 1 and 



The explicit forms of g' (ft) , g" (ft) are calculated as 



9' {ft) = 3 
9" (A) = 3 



C ° 8 Vcosh^V -IV-c-h 
^(^-l)ftanh^zA 



dy log (2 cosh-) 



tanh 



y/c + ft\ y/c 



(3.28) 
(3.29) 



S 



ceases to exist at |A| = 1. It is natural to ask what happens when |A| > 1. It is likely 
that at |A| = 1 the system undergoes the phase transition from conformally symmetric 
(deconfining) phase in |A| < 1 to broken (confining) one in |A| > 1. Let us briefly discuss 
the possibility of this phase transition by employing the results obtained above. 

For this purpose we consider the same system in the U(iV) (not SU(iV)) gauge group 
and take the diagonal U(l) holonomy in the thermal direction into account, which has the 
leading contribution in 1/N expansion. The self-energy and free energy density including 
the holonomy contribution can be respectively obtained directly from the self-energy ( |3.6| ) 
with ([3.15|) , (|3.14j) and the grand canonical potential density (|3.21|) by performing an 
analytic continuation for chemical potential \i — > ia, where a represents the diagonal U(l) 
holonomy. After this manipulation the solution of the self energy becomes 

f(p) = ^log (2(cosh ^p 2 + c + cosd)) , g(p) = - /(p) 2 , (4.1) 
e W = e ^ + + 2cosd, (4.2) 

where a — There exists a solution in < a < n when |A| < 1. In the region |A| > 1 
there also exists a solution in a c < at < n, where a c is a function of A taking values between 
it/ 2 and 27r/3.[] The free energy density is 



F{T,a) = 




Since the holonomy is dynamical, the large N leading free energy is determined by varying 
the holonomy satisfying (|4.2|) so that the free energy is minimized. It turns out by simple 
numerics that the value of the holonomy which minimizes the free energy is trivial when 
|A| < l|] but becomes non-trivial, a = a c , when |A| > 1. 

This result suggests that the first order phase transition may happen at |A| = 1. How- 
ever it also turns out that the minimized free energy (|4.3| ) at a = a c in |A| > 1 is positive, 
which means the pressure is negative and thus the system is no longer stable. This is 
presumably a signal of break-down of our assumption for the system to have conformal in- 
variance. Therefore operators neglected in this analysis, which break conformal symmetry, 
would play an important role to investigate the phase in |A| > 1. We leave detailed study 
in the future. 



5 Summary and discussion 

We have studied SU(iV) Chern-Simons theory coupled to a fermion in the fundamental 
representation with chemical potential for U(l) flavor symmetry and temperature in the 
't Hooft limit. We have calculated the self energy, the momentum distribution and several 
thermodynamic quantities exactly as a function of A, which suggests that the solution 

6 Precisely speaking, a c is given by a c (A) = arccos Max^ f| exp(-j^) — cosh yfcj , where M.asi x A{x) is 

the maximum value of A(x) with respect to x. 

7 This result was assumed in |J to determine the free energy. 
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exists up to |A| = 1 for a fixed value of chemical potential. We have discussed instability of 
the system beyond |A| = 1 by considering U(iV) gauge group and using the diagonal U(l) 
holonomy 

Even though the technique to solve Chern-Simons vector model seems to work well, 
there remains a puzzle against the result of calculation of three point functions flll2|-|l4| . The 
result of three point functions suggests that there exists a continuous parameter (just like 
the 't Hooft coupling) which interpolates between free/Gross- Neveu fermion vector models 
and critical/free scalar vector models respectively. However it seems this "bosonization" 
in three dimensions is not seen from the free energies computed by using the technique in 



H |19| . We do not know the answer of this puzzle. 

Related to this issue it is quite interesting to study the sub-leading correction in 1/N 
expansion since it might be possible that 1/N expansion breaks down in such a limit as 



pointed out in [14, 22 1. In that case the free energy obtained in 0, 19 would be modified 



so that non-planer diagrams are taken into account. At least we expect the next-leading 
contribution will be obtained by considering the contribution of the ladder diagrams, which 
will be technically much harder but considerably important to compute. 
We hope to come back to these issues in the future. 
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